Destroying superfluidity by rotating a Fermi gas at Unitarity 
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We study the effect of the rotation on a harmonically trapped Fermi gas at zero temperature under 
the assumption that vortices are not formed. We show that at unitarity the rotation produces a 
phase separation between a non rotating superfluid (S) core and a rigidly rotating normal (N) gas. 
The interface between the two phases is characterized by a density discontinuity tin/ws = 0.85, 
independent of the angular velocity. The depletion of the superfluid and the angular momentum 
of the rotating configuration are calculated as a function of the angular velocity. The conditions of 
stability are also discussed and the critical angular velocity for the onset of a spontaneous quadrupole 
deformation of the interface is evaluated. 

PACS numbers: 



The effect of the rotation on the behavior of a super- 
fluid is a longstanding subject of investigation in con- 
densed matter as well as in nuclear systems [l|, 0j ■ Due 
to the irrotationality constraint imposed by the existence 
of the order parameter a superfluid cannot rotate like a 
normal fluid. This implies the quenching of the moment 
of inertia at small angular velocities 51 and the occurrence 
of quantized vortices at higher fi. These peculiar features 
have been experimentally studied in a systematic way in 
superfluid helium and, more recently, in Bose-Einstein 
condensates (see e.g. Q). 

The recent realization of ultracold Fermi gases with 
tunable interaction close to a Feshbach resonance has 
opened new opportunities to investigate the effects of su- 
perfluidity along the so called BCS-BEC crossover (for 
recent reviews see, for example, 0, Q)- While on the 
BEC side of the resonance (small and positive values of 
the scattering length) the physics of a Fermi superfluid 
approaches the behavior of a Bose-Einstein condensed 
gas of dimers, on the opposite BCS side (small and nega- 
tive values of the scattering length) superfluidity is much 
more fragile because of the weakness of Cooper pairs. A 
further challenging regime is given by the unitary gas 
of infinite scattering length where the system exhibits a 
universal behavior and, in many respects, behaves like a 
strongly interacting fluid. Quantized vortices have been 
recently observed on both sides of the resonance _6] , in- 
cluding unitarity, providing a unifying picture of the be- 
havior of the system along the crossover. 

Previous theoretical work on superfluid rotating Fermi 
gases has mainh^ focused on the dynamics [3] and on the 
instability ^, ^ of configurations with high vortex den- 
sity. The purpose of this Letter is to investigate the be- 
havior of the gas when quantized vortices are not formed. 
Experimentally this scenario is achievable through an 
adiabatic ramping of the the angular velocity of the ro- 
tating trap starting from a configuration at rest. In the 
case of trapped BEC's this procedure has permitted to 
reach values of angular velocity significantly higher than 



the Fcynman critical angular velocity for the formation 
of quantized vortices whose nucleation is inhibited by the 
presence of a barrier. Under these conditions new phys- 
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FIG. 1: Upper panel: without rotation the system is fully 
superfluid. Lower panel: the rotation of the trap favours the 
formation of a phase separated state, with a superfluid core 
surrounded by a normal part exhibiting the bulge effect of the 
rotation {Q, — 0.45ui±). 

ical phenomena occur like, for example, the spontaneous 
breaking of rotational symmetry caused by the energetic 
instability of the surface modes [l^] and, at even higher 
angular velocities, the occurrence of dynamic instabili- 
ties jll| . These phenomena have been observed experi- 
mentally in BEC's ^12l] . confirming in a qualitative and 
quantitative way the correctness of the irrotational hy- 
drodynamic picture of rotating superfluids. 

The key point of the present work is that trapped su- 
perfluid Fermi gases can behave quite differently from 
BEC's because pairs are easily broken by the rotation. In 
the following we will restrict our discussion to the most 
relevant unitary gas. While at small angular velocities 
the superfluid is unaffected by the rotation of the trap 



|13| , we predict that at higher angular velocities the rota- 
tion results in a phase separation between a non rotating 
superfluid core and a rigidly rotating normal component 
(see Fig. (HI)). The mechanism of pair breaking is very 
intuitive. In fact near the border of the cloud, where the 
density n is small, the energy cost for destroying superflu- 
idity is also small, being proportional to n?'^. Viceversa 
the centrifugal energy gained by the normal phase can be 
large, being proportional to fl^R^ where R is the radius 
of the cloud. Notice that this mechanism cannot occur 
in the BEC regime due to the high energy cost needed 
to break the dimers. Furthermore, differently from the 
nucleation of vortices, the formation of the normal part 
is not inhibited by the presence of a barrier as all the rel- 
evant energy scales are vanishingly small near the border 
and the energy gain is ensured as soon as one starts ro- 
tating the trap. The appearance of a normal part due to 
rotation is not peculiar of ultra-cold atomic Fermi gases. 
For example, superfluidity is known to be unstable in 
nuclei rotating with high angular velocity [3|- Such an 
effect is also related to the so-called intermediate state in 
type-I superconductors, where the role of the rotation is 
played by the external magnetic field (see, e.g., [ij|). 

We are interested in a Fermi gas confined by a har- 
monic potential V{r) — m{ujlx'^ + uj^y^ + uj'^z^)/2, ro- 
tating with angular velocity il along z. We study the 
problem in the frame rotating with the trap where the 
potential is static and the Hamiltonian is H — flLz ■ In 
the local density approximation the energy of the rotat- 
ing configuration at zero temperature can be written as 



E= dv 



e{n) + V{r) + -mv'^ 



■mil{r X v)z — ^ 



(1) 

where e(n) is the energy density per particle, v the ve- 
locity field, n the chemical potential and n the density. 

By terming Rs{0, </>) the interface separating the super- 
fluid from the normal component, the integral in Eq.((T]) 
splits in two parts. The internal superfluid core occu- 
pies the region r < Rs{0,4i) and the surrounding nor- 
mal phase is confined to Rs{d, (j)) < r < R^id, 4>), where 
R^id, (f>) is the Thomas- Fermi radius of the gas where the 
density vanishes. The energy densities in the two phases 
are given, respectively, by 



es=es||^(6^^ns)^/^ 
5 Zm 



^ "-^ /c 2 n2/3 



(2) 



(3) 



where ng (n-^) is the superfluid (normal) density and 
the dimensionless parameters ^g — 0.42 and ^n = 0.56 
account for the role of interactions in the two phases. 
Their value has been calculated in [1^, ^^ employing 
Quantum Monte Carlo simulations. 

The equilibrium is found by minimizing the energy 
with respect to the densities and the velocity fields of 



the superfluid and normal part, as well as with respect 
to the position of the border surface. Notice that this pic- 
ture ignores surface tension effects, a plausible assump- 
tion in the limit of large samples. The superfluid velocity 
obeys the irrotationality constraint and can be written as 
vs — V<I>. Varying the energy with respect to the veloc- 
ity potential $ yields the continuity equation 



V- ((V$- J7x r)ns) = 0, 



(4) 



while variation with respect to the superfluid density ng 
leads to 



M = &^(67r2«s)2/^ + yS(r), 



(5) 



where F^(r) — V{r) + ^niv'^ — mn{rx-vs)z is the effective 
harmonic potential felt by the superfluid. 

Using the same procedure for the normal part without 
the irrotationality constraint we have vn — ft x r, i.e. 
the normal component rotates rigidly. The variation with 
respect to the normal density yields an equation similar 
to dSl) 



fi2 



M = eN:^(6^^nN)'/' + W"(r), 
Zm 



(6) 



where the effective harmonic potential y'^(r) is now 
quenched by the rigid rotation according to (i-i^^)^ = 



N\2 



n' 



By varying the energy ([1]) with respect to RsiO, 0) we 
eventually find the equilibrium condition for the coexis- 
tence of the two phases in the trap. The resulting equa- 
tion implies that the pressure of the two phases be the 
same: n'^{des/dns) = n'^{de-f^/dn-f^). Using Eqs.© and 
([5|, one then predicts a density discontinuity at the in- 
terface given by 
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0.85, 



(7) 



independent of the angular velocity. The equal pressure 
condition ([7]) combined with Eqs.([5]) and ([6]), results in 
the useful relationship 



i^,-V'{r))^Ji^,-V''ir)), 



(8) 



which determines the surface RsiO, (p) separating the su- 
perfluid and the normal part. 

In the following we assume uix = i^y ^ loi_ and we 
consider the solution vs = 0, corresponding to a non 
rotating axi-symmetric superfluid and hence to V^ — V. 
In this case wc find 



RU9) ^^{cicos^ 
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(9) 

where 6 is the polar angle. On the other side the Thomas- 
Fermi radius i?N of the normal gas is fixed by the condi- 
tion /i — V yielding 

R^{9) = ^ {ujI cos^ e + Lul sin^ 9 - fl^ sin^ ^T^ > ^s(^) 

(10) 



The value of fi is fixed by the normaUzation condition 



rigdr 



r<Rs{e) 



nyidr = N 



(11) 



Rs(0)<r<Rt^{9) 



where N is the total number of particles. While for 17 = 
the system is completely superfluid, for fi > it phase 
separates into a superfluid and a normal component char- 
acterized by the density jump ([7]) at the interface. This 
behavior shares interesting analogies with the the phase 
separation between a superfluid and a normal component 
exhibited by polarized Fermi gases [13, \lB, U^ where a 
jump in the density at the interface is also predicted to 
occur [20| ■ By tomographic techniques [l^ it is nowadays 
possible to measure the density "in situ", thus the pre- 
dicted discontinuity, and hence the value of ^s/^n, should 
be observable experimentally. 

The radii Rs and _Rn coincide at 6* = which means 
the absence of the normal part along the z-axis. In the 
plane of rotation {6 — tt/2) the difference between the 
two radii is instead maximum and becomes larger and 
larger as fl increases. In particular the radius of the 
superfluid is always smaller than the radius of the cloud 
in the absence of rotation, while the radius of the normal 
gas is always larger due to the bulge effect produced by 
the rotation. 
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FIG. 2; Density profile for 9 — ■k/2 of the rotating Fermi 
gas at fi = 0.4:5uj± (full line). The profile in the absence of 
rotation is also shown (dashed line). Inset: Superfluid radius 
versus Q/lu± . 

In Fig. ([2]) we plot the densities Tig and tin as a function 
of the radial coordinate at 6 — tt/2 in a spherical trap 
for fi = 0.45w_L. The densities and the radial coordi- 
nate have been renormalized with respect to the central 
density rig and the Thomas-Fermi radius Rq of the su- 
perfluid at rest. The inset shows the superfluid radius Rs 
renormalized by Rq as a function of the angular velocity 
Q/u)±. 

From the knowledge of the density proflles and from 
the radii Eqs.© and (fTU]) we can calculate the number of 
particles in each phase. In Fig. ([3]) we show the ratio be- 
tween the number of particles iVg in the superfluid phase 



and the total number N as function of the angular veloc- 
ity. The higher the angular velocity, the more particles 
prefer to stay in the normal phase and thus the superfluid 
is depleted. At small angular velocities the depletion of 
the superfluid follows the law Ns/N = 1 - (j^)^^^ n^. 
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FIG. 3: Depletion of the superfluid as a function of the trap 
angular velocity for an axi-symmetric configuration (vs. 

Another important observable is the angular momen- 
tum Lz- For an axi-symmetric conflguration the su- 
perfluid does not carry angular momentum which is 
then provided only by the normal component: Lz = 
il J dr {x'^ -|-y^)nN. The total angular momentum then 
increases with fl and eventually reaches the rigid body 
value at 11 = uj± (see Fig.((31)). The angular momentum 
of a rotating configuration has been measured in EEC's 
by studying the precession phenomena exhibited by the 
surface excitations j21| . 




0.2 0.4 0.6 0.8 



Q/O)^ 



FIG. 4: Angular momentum in units of the rigid value as 
a function of the trap angular velocity for an axi-symmetric 
configuration. 

While for small values of the angular velocity {fl < 
0.2uj±^) the superfluid is robust, it is remarkable that even 
at angular velocities far from the centrifugal limit the de- 
pletion of the superfluid and hence the angular momen- 
tum of the system are sizable. 



A major issue concerns the conditions of stability of 
the rotating configuration discussed in this work. Let us 
first consider the question of energetic stabihty. We have 
shown that, in the frame rotating with the trap, the phase 
separated configuration is energeticaUy favoured with re- 
spect to the configuration where the whole gas is at rest 
and superfiuid. This is true for any value of the angu- 
lar velocity. When il exceeds a critical value of order 
{h/mB?) ln{R/d), where d is the healing length fixed, at 
unitarity by the interparticle distance, quantized vortices 
become an even more favourable configuration. This en- 
ergetic instability is not however expected to be a severe 
difficulty if one increases the angular velocity in an adi- 
abatic way because the presence of a barrier inhibits the 
access to the vortical configuration as proven experimen- 
tally in the case of EEC's [22]. At higher angular veloc- 
ities the axi-symmetric superfiuid configuration is even- 
tually expected to exhibit a surface energetic instability, 
undergoing a continuous shape deformation, similarly to 
what happens in Bose-Einstein condensates [13, 113, ^M ■ 
This effect is accounted for by the solution vg ^ of 
Eq.Q, corresponding to a spontaneous breaking of ro- 
tational symmetry in the superfiuid component. In the 
case of Bose-Einstein condensates a quadrupole instabil- 
ity takes place at flcr = iLj±/V2 and for larger values of ft 
the solution vg = corresponds to the so called overcrit- 
ical branch J^] . A different value of flcr is predicted for 
the rotating Fermi gas, due to the new boundary condi- 
tion imposed by the presence of the normal component. 
This condition requires that the velocity of the super- 
fluid be tangential to the interface. In the simplest case 
of a 2D configuration we find the value ficr = 0.45w^ for 
the emergence of a spontaneous quadrupole deformation. 
The experimental measurement of flcr would provide a 
crucial test of the consequences of the phase separation 
caused by the rotation of the unitary Fermi gas. 

An even more challenging question concerns the emer- 



gence of dynamic instabilities. In the case of a rotating 
BEC a dynamic instability takes place at values of il 
slightly larger than uj±/^/2 and corresponds to the ap- 
pearance of imaginary components in the frequency of 
some hydrodynamic modes [ll| . In the case of the rotat- 
ing Fermi gas discussed in this work a dynamic instability 
might be associated with the Kelvin-Helmholtz instabil- 
ity of the interface between two fiuids in relative motion 
(see.e.g [2J]). However, if the densities of the two flu- 
ids are different, an external force stabilizes the two-fluid 
system against the appearance of complex frequencies in 



the low energy excitations of the interface [2j| . This is 
actually our case where the density of the two phases 
exhibits the gap ([7]) and the system feels the external 
force produced by the harmonic confinement. We conse- 
quently expect that the system be dynamically stable at 
least for moderately small values of the angular velocity. 

In conclusion we have shown that an ultracold Fermi 
gas at unitarity can separate into a superfiuid and a nor- 
mal component as a consequence of the adiabatic ramp- 
ing of the trap rotation. The formation of the rotating 
normal component requires that the trap transfers angu- 
lar momentum to the gas within experimentally accessi- 
ble times. Its realization would open the unique possibil- 
ity of exploring the Fermi liquid behaviour of a strongly 
interacting gas at zero temperature. Important effects 
to investigate are, for example, the zero sound nature of 
the collective oscillations and the behaviour of viscosity. 
The detailed study of the energetic instability associated 
with the spontaneous breaking of rotational symmetry as 
well as the effects of the rotation on the polarized phase 
{N-^ ^ Ni) will be the object of a future work. 
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